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, ^ Abstract 

Active Brownian Particles are self-propelled particles that move in a dissipative 
medium subject to random forces, or "noise". Additionally, they can be confined by 
an external field and/or they can interact with one another. The external field may 
actually be an attractive marker, for example a light field (as in the experiment) or 
an energy potential or a chemical gradient (as in the theory). The potential energy 
can also be the result of interparticle attractive and/or repulsive forces summed 
over all particles (a mean field potential). Four, qualitatively different motions of 
the particles are possible: at small particle density their motions are approximately 
independent of one another subject only to the external field and the noise, which 
results in moving randomly through or performing rotational motions about a cen- 
tral point in space. At increasing densities interactions play an important role and 
individuals form a swarm performing several types of self-organized collective mo- 
tion. We apply this model for the description of zooplankton Daphnia swarms. In 
the case of the zooplankton Daphnia (and probably many other aquatic animals 
that form similar motions as well) this vortex is hydro dynamical but motivated 
by the self-propelled motion of the individuals. Similar vortex-type motions have 
been observed for other creatures ranging in size from bacteria to flocks of birds 
and schools of fish. However, our experiment with Daphnia is unique in that all 
four motions can be observed in controlled laboratory conditions with the same an- 
imal. Moreover, the theory, presented in both continuous differential equation and 
random walk forms, offers a quantitative, physically based explanation of the four 
motions. 

Key words: Active motion, self-propelled particles, vortex motion, noise-induced 
bifurcation, symmetry breaking, Zooplankton 
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1 Introduction 



The theory of Active Brownian Particles (ABP) has been developed in some 
detai l over the past several years and published largely in physics litera- 
ture (Schimanskv-Geier et all Il99^ : lErdmannI. Il997l: [Schweitzer et al. . Kl 



Eheling et alllmwAlKnlinanii el nlll'jOOdHScliweil /eil.b()(d). The ABP theory 



is based on continuous motions in time and is described by a set of stochastic 
differential equations. A second approach has been developed, which is based 
on the random walk theory (called RWT from now on). It has been applied 
to p articles that moye in di screte "hops" with a pause at the end of each 
hop ( Ordemann et ah L 2003alfb c ) . Though or i ginall y conceived without a self- 



propelling velocity ( Schimanskv-Geier et al. . 1995j ). the present ABP theory 



deals with particles that are self-propelled and move through a dissipative 
medium, such as water. The particles expend energy in order to move and to 
supply their own metabolic requirements, while foraging for food represents 
an energy uptake. 

These notions are conveniently expressed in terms of an effective, velocity de- 
pendent dissipation, 7(f), which can assume positive and for small velocities 
possibly negative values. F or 7(f) < 0, energy i s pumped into the particle and 
may be stored in a depot (|Ebeling et all Il999l ) (there is more than sufficient 
food to supply the energies of motion and metabolism), while for ^(v) > the 
rate of energy uptake and conversion is insufficient to maintain the supply in 
the depot causing the particle to slow down. The particle then finds a velocity 
for which, on average, the effective dissipation is zero, and this value is the 
self-propelling velocity. See Fig. 1 below. The right/left symmetry of these 
motions may be broken by parallelizing interactions, what could be any kind 
of a velocity aligning force. Single or groups of ABPs can move in a central 
potential which insures that they are confined to some region in space. Here 
we confine our model to a quadratic potential representing a linear restoring 
force. In addition, a random force, or noise, is applied in order to represent the 
variability inherent in a real animal's motions. In contrast to the ABP theory, 
the RWT begins with the assumpti on of an average, c onstant self-propelling 
velocity, as does another approach ( Vicsek et al. . 19951 ). Each particle moves 



by hopping and turning. At the end of each hop a change in direction is taken 
from an experimentally measured turning angle distribution and a "kick" to- 
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ward a center of motion is applied. Random selections from the turning angle 
distribution represent the noise and the kick insures that the particles are con- 
fined. Here we take the kick to be a force linearly proportional to the distance 
between the particle and the center of the field. Both theories describe systems 
that operate far from thermodynamic equilibrium, since motion is maintained 
in the dissipative medium by a continual uptake of energy. 



In the ABP theory, coherent motions at low density are modeled by interac- 
tions with the central field plus noise only. Two attractors motions are sep- 
arated b y a bifurcation control led by the magnitude of the self-propelling 
velocity ( Schweitzer et allll998f ): one a noisy fixed point and the other a pair 
of symmetric, limit cycles (recurrent closed paths) with equal probabilities to 
show a clockwise (CW) or counterclockwise (CCW) rotational sense in a plane 
fjErdmann et all l2000li . Se e Fig. 2 below Other systems, for example molecu- 
lar motors (|Juficher et all Il997l ). which are also far from equilibrium and be- 
have according to a similar velocity dependent dissipation, also show bimodal 
velocity d istributions, that is, m otion in two opposing directions that is equally 
probable ( Badoual et al. . 2002). The theory presented here is two-dimensional 
so that these motions are confi ned to the x-y plane. A particle-particle attrac- 



so mat tnese motions are conn ned to tne x-y plane. A particle-particle attrao 
five in teraction is then built in ( Ebeling and Schweitzer! . 200lt Erdmann et al. 



2002a)). The interaction is global in the sense that every particle feels the at- 



traction of every other one but the strength decreases with pair separation 
distance. Such interactions result in a central mean field and predict swarm- 
ing behavior, the gradual increase of particle density around the symmetry axis 
of the central field. Finally, a particle-particle hydrodynamic interaction is in- 
cluded that models the tendency for nearby particle s (feeling one another's 
flow field) to align the directions of their velocities (jErdmann and Ebeling . 
20031) . This interaction breaks the symmetry of the two limit cycles leading 



to a collective vortex-type motion in one or the other rotational direction. 
A short range repulsive potential that results in avoidance when two parti- 
cles come closer than some minimum distance can also be used to break the 
symmetry and lead to vortex formation. The sense of rotation of the vortex 
depends on an initial condition that is random, so that the final motion can 
be either CCW or CW. 



Thus the ABP theory predicts four basic types of motion if the particles are 
confined: 1) noisy fixed point whereby a particle moves randomly about and 
through a single point in space; 2) two limit cycles whereby the particle moves 
on closed paths with CCW or CW directions equally probable. The second 
type can be separated again in two different types of motion if we investigate 
larger populations of interacting individuals: 3.1) swarming whereby initially 
randomly and widely dispersed particles first move toward a spatial center and 
then start rotating around this center (rotational motions in both directions 
are equally probable, and thus the net rotational motion is on average zero); 
and 3.2) vortex whereby all particles spontaneously align their velocities and 
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commence rotation in the same direction, because of a higher density of the 
individuals. As we show below, all four of these motions can be observed in 
laboratory experiments within a single species: the light sensitive zooplankton 
Daphnia. 



Such motions, in addition to other collect ive behaviors, ar e commonly o b 
served dParrish and Edelstein-Keshet , 1999h in schools of fish (IWeihsLll973allbl: 
Hall et al. . Il986l : IParrish et all 20021) . bacterial colonies (iBen- Jacob et al 



and flocks of birds ( Larkin and FraseL 1988| ) . Vortex motions 



1994), slime molds (Ibevine and Revnoldlh 99ltlBorrneTJ . il 99SllDao et alH 2000) 



were s 



cally observed long ago in colonies of an ts by Schneirla (lAronson et al 



ecifi- 



Il999l) 

Okubo and Levin .12001 : IParrish et al. , 20021 ) . The idea of a self-propelling ve- 



and more recently in bacterial colonies (|Czir6k et all 119961 ) and slime molds 
( 



Rappel et all I1999T) . Swarming is of course q uite well known. See for example 



locitv arose early ([Weihsl . 1973b[ ). and was expressed by Sakai whose model pre- 



dicted the aforementioned four types of mo tion already in 1973. See a discus- 
sion on that early work in the useful book ((Okubo and Levinl . 120011 1 . General 



reviews on the physics of individual part i cle mo tions and aggregation ( swarm- 
ing) has been given bv lBen-.Tacob et al.1 (|2000f ) and iFlierl et all f)l999l ). More 



recent theories, based on the idea of a self-propelling velocity, incorporated 
various coupli ngs among the particles such as directional averaging of the vec 



tor velocities (IBen- Jacob et al. . 1994 : Vicsek et al. . 19951 : Levine et al. . 2001 



Couzin et alll2002|) . A hydr odynamic generalization of t hese "rule based" the- 



ories was later developed bv lToner and Thl (|l995l ll998). In contrast to these 



the ABP theory, based on the energy depot model ([Ebeling et al 
represented by stochastic differential equations. 



1993), 



is 



Our objective here is to summarize the ABP and RWT theories without bur- 
dening the reader with excessive detail and to compare the predictions of those 
theories at each step in their development with our experimental results on 
the motions of Daphnia. 



2 Outline of the ABP Theory: single particle case 



In 1827 the British botanist, Robert Brown, observed irregular motions of 
microscopic particles in water. He noted that their motions resembled that of 
living creatures. Later all such entities whose motions are mediated by random 
forces came to be called "Brownian Particles". In a two-dimensional spatial 
plane, x-y, they are located by a vector r, whose magnitude is r = ^/x 2 + y 2 
and move with velocity, d t f = v. Physical models represented by stochastic 
differential equations, called Langevin equations, describe the motion of these 
particles. They represent a balance of forces on the particle: the random force, 
or noise V2D£(t), where D is the noise strength; the confining restoring force, 
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VL r (r), where U(r) is the external potential energy, the passive friction — 70^ 
and, in our case, the self-propelling force de(t)v that results from conversion 
of energy out of an internal energy depot e(t) into motion; 

md t v = -7o# + d e{t)v - VU(r) + y/2D£{t) . (1) 



The noise, £(t), in this equation is external in the sense that it is contributed by 
the environment in which the particle finds itself. In the case of ABPs this noise 
could arise, for example, from thermal collisions, hydrodynamic fluctuations 
or turbulence. This noise is assumed to be "white" , that is its autocorrelation 
function is a delta function, = 2DSijS(t — s), Gaussian distributed 

and with zero mean. As an identifier, we will call this the "external" noise. 
Further, from now on we will choose units in which m = 1. The time evolution 
of energy stored in the depot is given by e(t), which is described by a second 
equation, 

d t e(t) = q -ce(t)-dv 2 e(t), (2) 



that represents a balance among three processes: go which is the energy in- 
flow available to the depot, here called the uptake, for example from for- 
aging for food; —ce(t), which is a constant drain representing, for example, 
a metabolic requirement; and dv 2 e(t), which is the power dissipated neces- 
sary to propel the particle through the medium. Here, c and d are constants. 
We take the available energy inflow qo, to be uniformly distributed over the 
two-dimensional space, though in more general treatments it can be irregu- 
larly distributed, q(r) representing, for example , food availability in patches 



( Steuernagel et al. . 1994 : Schweitzer et al. . 1998| ) 



A key issue is the velocity dependent dissipation, 7(f). Here we mention two 
different commonly used treatments. The first is the Rayleigh law, 7(f) = 
—71 +72^ 2 , where the two 7's are constants serving to delineate constant and 
velocity dependent contributions to the dissipation. After initial transients 



die out, the particle assumes a self-propelling velocity, 



71/72. This law 



was originally d eveloped to des cribe dissipative energy pumping in systems 
involving sound ((Ravleighl . Il894r ) an d has been used to describe self-propelling 
velocities ap plied to fish schooling (iNiwal 11994 [ 19961 ). The second is a law 
developed bv lSchienbein and Grulerl (J1993), 7(f) = 70(1 —Vq/v), developed to 
describe the motions of va rious motile biological cells . The same friction law 
was used independently by Helbing and Molnarl (|l995| ) to describe the motion 
of pedestrians. Both of these laws have singularities, the former in a range 
of v ^> vq, and the later for v — > 0. A key simplification can be obtained by 
assuming that the energy uptake rate e{t) is slow compared with the dissipative 
rates. Equation (2) can then be adiabatically eliminated, dte(t) = 0, which, 
together with Eq. (1), results in a third dissipation law that does not suffer 



5 



from singularities and thus is more physically reasonable (lErdmann et al 
200nh . 



j(v) 



7o - 



god 

c + dv 2 



(3) 



where 70 is the dissipation constant in the limit of high particle velocity. This 
function is shown in Fig. 1 where the velocity ranges of energy pumping into 
the depot (super-critical pumping) and dissipation from it (sub-critical pump- 
ing) are indicated separated by the self-propelling velocity at zero effective 
dissipation, 



To 



c 
d 



(4) 



y 




-0.5 



Fig. 1. The velocity dependent friction function, Eq. (3) showing sub- and su- 
per-critical energy pumping regions separated by the self-propelling velocity for 
which the effective dissipation is zero. The Rayleigh and Gruler formulae are shown 
by the dotted and dashed lines, respectively. The values of the parameters are, 
q = 10, c = 1.0, 70 = 2.0, and d = 0.7 

The solutions of Eq. (1) using the formula Eq. (3) give the first two types of 
motion separated by a bifurcation. The bifurcation parameter can be defined 
as the ratio of constants, 



%d 

C70 



(5) 
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Introducing a force VU(f) = ou^f, attraction by a light shaft, chemotaxis 
and similar causes can be modeled. For such a quadratic attractive potential 
(the force is linear in r) two-dimensional motion with natural frequency ujq 
on closed paths in the x-y plane is possible. For /i < 1, there is sub-critical 
pumping of energy and the motion is random around the symmetry point 
of the central potential, called here the "noisy fixed point"; and for /i > 1 
the pumping is super-critical, and the resulting motion is a symmetric pair of 
noisy limit cycles. Particles move on the limit cycles with the self-propelling 
velocity \v\ > 0. These limit cycles are closed curves in the x-y plane, and 
their symmetry represents the fact that the two senses of rotation, CCW and 
CW, occur with equal probability. This bifurcation is shown in Fig. 2, where 
the probability densities of the velocity components v x and v y are plotted 
showing the monomodal peak (noisy fixed point) in Fig. 2(a) for /i < 1 and 
the bimodal density (limit cycles) in Fig. 2(b) for /i > 1. The insets show 
example trajectories in the x-y plane. 




(a) fi = 0.35 (b) n = 3.5 



Fig. 2. The probability densities of velocity in the two-dimensional plane showing 
the noisy fixed point for fi = 0.35 in (a) and the limit cycles for fi = 3.5 in (b). 
The other parameters are 70 = 2.0, D = 2.0, c = 1.0 and qo = 10. The insets show- 
typical trajectories in the plane 



The foregoing theory applies to a single particle or to a collection of iden- 
tical non-interacting particles moving at a unique self-propelling velocity in 
a attractive central potential. Below, however, we compare the theoretically 
predicted motions to those we actually measure for populations of Daphnia. In 
this case, the individuals are not identical. One way this population variability 
could be modeled is by building in a distribution of power conversion rates 
as given by the last term on the right in Eq. (2). In our populations, some 
Daphnia are young and presumably cannot convert energy at the same rate 
and thus cannot swim as fast as older individuals. We thus modify Eq. (3) by 
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allowing the strength of the power dissipated in motion to become noisy, 



d — > d + r)(t) 



(6) 



where r](t) is noise, which, in contrast to the external fluctuation, £(£), is 
an internal fluctuating variable within the population, meaning that different 
individuals have different values of d. We call this the population noise. It 
is Gaussian distributed with zero mean , has standard deviation a^, and is 
delta correlated (see also lErdmarm et al.l M , where a factor describing the 
efficiency of conversion is added to the model). This translates into a noisy 
dissipation as given by Eq. (3). Using Eqs. (l)and (3) with (6), we can obtain 
trajectories and measures directly comparable to the experimentally measured 
ones as we show in the next section below. Here we calculate the probability 
densities of velocity for two conditions of noise as shown in Fig. 3. The solid 
line represents density with the noise in Eq. (1) only (the external noise) and 
the dashed line represents both noises accounting for external and population 
variability. The net result of including population variability in the theory is 
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Fig. 3. The probability density of radii and velocities projected on the x-axis for 
the single noise £(t) in Eqs. (1) and (2) shown by the solid line, and for both noises, 
and r](t), in Eqs. (1) and (6) shown by the dashed curve. 

that there is now a distribution of both radii and velocity at which the particles 
travel, a quantity that can be directly measured in the experiment. Note that 
for some of the creatures the energy conversion is not efficient enough to reach 
the limit cycle. These participate in Brownian motion around the fixpoint 
(see inset of Fig. 2(a)). This is indicated by the maximum around zero in both 
graphs of Fig. 3. 



3 About Daphnia 



Daphnia are found in many species in fresh water over the entire globe. Fig- 
ure 4 shows a single individual. They have a long evolutionary history and 
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Fig. 4. A single Daphnia magna with a store of eggs under the carapace. The two 
large antennae enable the animal to swim. An array of smaller appendages enables 
feeding. The animals range in size from about 3 to 5 mm in length. 

consequently have developed a range of complex behaviors that insure their 
survival in various environments both stressed and nonstressed. Daphnia are 
known to heavily depend on phototaxis, for example whe n searching for food 



( Young and Gettv , 1987 ) , to confuse an d avoid predators ( Jensen et al. L 19981) 



or for performing diel vertical migration (|Zaret and SufferrJ . ll976tlRhode et al. 
2001), a behavior whereby Daphnia swim near the water surface while feeding 
at night, but return to the depths during the day in order to avoid visual 
predators (most planktivore fish). They are strongly attracted to light in the 
visible range (VIS), repelled by ultraviolet (UV) and blind to infrared (IR), 
three facts that we will make use of in the experimental design outlined below. 



It is n ot likely that Daphnia c an for m an image with their eyes accord- 
ing to iBuchanan and Goldbergl ((1981). Unlike bir ds and fish, no direct v i- 
sual alignment between Daphnia has been detected (jOkubo and Levinl . l200lh . 
They can, however, determ ine wavelength, intensity, an d direction of light 

( Smith and Macagnol 1990f ). Daphnia also use chemotaxis ( Larsson and Dodson . 

1993) and can detect water motions with their mechanoreceptors (|Haurv and Yamazaki . 
1995). Because the re is no known mechanism of long-range communication 
between Daphnia ( Larsson and Dodsorl Il993h in order to display coherent 
motions t hey must be individ ually attracted to a landmark such as a shaft of 
VIS light ( Jensen et al. . 1998| ). This can be easily demonstrated with a vertical 
flashlight beam in a dark aquarium as shown in Fig. 5. The swimming behavior 
of Daphnia is dominat ed by the fa ct that they live in a relatively low Reynolds 
number environment ( Zaret . 1980i ). They move with an average 'hopping' rate 
of approximately three moves per second and their overall sp eed is 4-16 mm/s 
with a sinking rate of approximately 3 mm/s ( Dodson . 1996f ). See an example 
3-dimensional trajectory in Fig. 6. Daphnia at low density (about one individ- 
ual per liter), swimming in the dark, do not show evidence of any long-range 
interactions. (There is, however, a short-range repulsion leading to avoidance 
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Fig. 5. Daphnia are attracted to VIS light as shown by this vertical flashlight beam. 
Three side view snapshots taken at 20 second intervals show the gathering of Daph- 
nia (white dots) in the shaft of light. 



Fig. 6. Typical trajectory measured on a Daphnia in darkness in three dimensions. 
This trajectory illustrates the hop and sink motion with a change of direction be- 
tween successive hops. 

maneuvers during close encounters.) Therefore, their motion can be modeled 
by the ABP theory of single self-propelled particles moving in a central field 
and subject to noise. 



4 Experiments with small Daphnia density 

The apparatus designed to investigate the aforementioned motions is shown 
in Fig. 7. It consists of a rectangular aquarium (50 x 20cm, water level 25cm) 
with a vertical translucent plastic tube (30 cm long by 1 cm diameter) mounted 
in the center. The tube is illuminated with VIS light from above by a fiber optic 
cable thus forming a shaft of light propagating radially outward. This forms a 
vertical marker with cylindrical symmetry to which the Daphnia are attracted. 
A filter to remove IR from the light shaft is placed between the illuminator 
and the fiber optic cable. Daphnia within the aquarium are illuminated by 
IR from arrays of filtered IR-Light Emitting Diodes located on either side as 
shown. IR sensitive video cameras (Baxall, Cohu) fronted by filters to remove 



x 



z 




y 



10 



Cameras. 
V!S-biock 
Filters „ 



D 




VfS Ught 
source 

!R 

tHuminators 




■ Aquarium & 
Daphnia 

mirror 




Fig. 7. The apparatus, showing the light source and cylindrical VIS light shaft, 
the IR illuminators and cameras. Daphnia motions can be observed and recorded 
simultaneously in side and bottom views. Appropriate filters are arranged such that 
Daphnia are recorded in IR while responding to and seeing only VIS. 

VIS light obtain side and bottom views of the aquarium. The bottom view is 
afforded by a 45° first surface mirror. Images from the cameras are combined 
on a split screen viewer, recorded simultaneously on magnetic tapes, then 
later digitized and saved on large capacity disks in a PC. The trajectories are 
analyzed with the tracking software (Chromotrack@ Version 4.02 from San 
Diego Instruments, or Tracklt from IguanaGurus) . In a later improvement, the 
video cameras were replaced by an IR capable camcorder (Sony DCRTRV80) 
with niters to block VIS. With this arrangement digital data can be streamed 
directly to the computer. 



After placing 20 to 80 Daphnia in the water with algae (Scenedesmus quadric- 
uld) for food and leaving them to acclimate and to distribute uniformly for at 
least 15min, the VIS light shaft is switched on. Being individually attracted 
to the optical marker, a substantial number can be observed to move on cir- 
cular paths around the light shaft in both directions (CW and CCW) with 
frequent reversals of direction. The radii of their tracks are large enough (typ- 
ically 1 to 5 cm) to exclude the possibility that this behavior occurs simply 
due to the hydrodynamic sensing of an artificial object presumably with their 
mechanoreceptors. Experiments show that Daphnia can sense static objects at 
dista nces not exceeding about one to two body lengths (|Haurv and Yamazaki . 
1995j b Similar circling behavior was observed in the absence of any solid object 
in more recent experiments using a high intensity flashlight shining vertically 
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into the water instead of the solid light tube ( Nihongi et all l2003f ). Typical 
tracks are shown in Fig. 8(a) for about 40 individuals and Fig. 8(b) for a single 
individual that makes seven complete CCW turns and a reversal of direction 
around the central marker. Arrows indicate the rotational directions. The tra- 
jectories look as in the theory shown in the inset of Fig. 2(b). While these 
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Fig. 8. Bottom view of individual Daphnia motions around the light shaft, (a) About 
40 distinguishable individuals, visible as a chain of five black dots for each animal, 
(grey scale inverted for easier recognition, time interval between single dots 0.3 sec) 
are engaged in circular motions indicated by filled arrows with approximately equal 
numbers traveling CCW and CW. The 1 cm diam. white circle is the end view of 
the light shaft with VIS light blocked. The dashed arrows in the outer areas of the 
field of view indicate Daphnia that are moving towards the light, (b) The track of 
a single individual making seven CCW turns (solid line, start marked by arrow) 
around the light shaft and showing a reversal toward the end of the track (dashed 
line, end marked by arrow, coordinates in mm). 

results from experiment and theory seem quite similar, statistical measures 
are necessary for quantitative comparisons. 

An essential point is that single Daphnia make rotational motions around the 
central attractant. This indicates that, for Daphnia at least, circling is not an 
emergent or self-or ganized motion that occurs a s an in herent property of a 



swarm of animals (jParrish and Edelstein-Keshetl . Il999|h Circular motions in 
zooplankton at low densi ty were observed previously ((Young and Gettyl . Il987 ; 
Young and Tavloi . 1990| ) but only incidentally in the course of different ex- 



periments. Below we show that high-density swarms of Dayhnia are c apabl e 



of vortex motions similar to those observed in bacteria (Czirok et all 1199 



slime molds ( Rappel et al. . 1999^ and fish ( Parrish et al. . 2002| ). But in Daph- 



nia the onset of the vortex is the result of a symmetry-breaking process that 
transforms the two symmetric limit cycle motions (CCW and CW), observ- 
able for single individuals in both theory and experiment, into motion in one 
direction only. Certainly individual birds in some cases follow rotational paths 
around central attractive markers ( Larkin and Frasei Il988h . Possibly bacte- 
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ria and fish do also. It is therefore not demonstrated that the vortex motion 
sometimes observed in swarms of these animals is an "emergent" property 
unrelated to the motions of individuals. 



We turn now to statistical measures of Daphnia motions at small density. The 
dominant part of the motion of single Daphnia in the light field takes place 
in the horizontal plane. Indeed individuals tend to circle in layers with little 
vertical migration, as we have o bserved in 3-dimen sional tracks obtained in 
the laboratory of J. R. Strickler ( Nihongi. Ai . 2002[ ). It is therefore sufficient 
to restrict the characterization of the observed circling behavior to two dimen- 
sions, analyzing only the bottom-view pictures of the aquarium. To quantify 
the amount of circling, we measured the probability distribution P(6) of the 
heading angles for several Daphnia observed circling around the light shaft. 
The heading angle, 9, is defined as the angle between the vector direction of a 
single hop and the direction to the light shaft measured in the middle of each 
hop as shown in Fig. 9(a). This distribution is double peaked with maxima lo- 



light shaft 




0.4 

^-s 0.3 

0.2 

0.1 



-160 -80 



80 160 



6 
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(b) 



Fig. 9. The heading angle, (a) Definition of the heading angle, 8, as measured 
for every hop movement of an individual, (b) Histogram of heading angles, P(0), 
obtained from 624 hops f rom the tracks of four individuals. Plot (b) is adapted from 
Ordemarm et~aTI {2n03hl ). 

cated at approximately ±90°. It indicates two properties of the motion: First, 
the twin peaks at ±90° demonstrate that the average motion is resolved into 
two rotational directions, CW (+90°) and CCW (—90°). Second, the magni- 
tudes of the two peaks are nearly equal, demonstrating that the CW (+90°) 
and CCW (—90°) motions occur with approximately equal probabilities. Both 
of these experimental observations are predicted by the ABP theory. 

For the same set of hops we obtained a histogram representing the probability 
distribution P(r) of the distance r of the Daphnia to the light shaft as shown 
in Fig. 10(a). Figure 10(b) shows a histogram of the angular momentum, again 
with twin peaks of approximately equal magnitudes. In Fig. 10(c) we show a 
histogram of the average length of circling in one direction (CW or CCW) 
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before reversing the direction. Further on we determined the average length 
of ' 
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Fig. 10. Histograms of (a) the circling radii, P(r), (b) the angular momentum, 
P(L ang ), and (c) the directional reversals, P{Mq). These data were obtained from 
the same set of hops usin g the same track data as for Fig. 9. Plots are adapted from 
lOrdemann et al.l feOQShl lT 



5 The Many Particle ABP Theory 



When there are many particles with a long range attractive interaction be- 
tween them, the global effect is that the particle-particle interactions can be 
replaced by a mean field. This field is the result of a central potential that 
causes confinement of the particles. This phenomenon can actually be demon- 
strated in experiments with the zooplankton Chydoridae. A long range attrac- 
tive interaction can be switched on by causing one or a very few individuals to 
become visible to all the others by shining light on them. When that happens 
all animals are attracted to the visible few. The light scattered from the few 
must, however, dominate the visual environment. This demonstration has been 
realized in the laboratory of J. R. Strickler by directing a narrow, blue, VIS, 
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laser beam vertically into an aquarium which, apart from the beam, is dark. 
The beam itself cannot be seen by the animals in the tank since the medium 
is clear and the beam is very narrow. Many animals are initially distributed 
either at the bottom or throughout the aquarium. One or a few individuals by 
accident swim into the beam, become brightly illuminated, and thus become 
visible to all the rest. The remaining individuals are attracted by this light, 
immediately begin migrating toward the beam and become illuminated them- 
selves, thereby intensifying the attraction, a positive feedback effe ct whose 
streng th is dependent on the number illuminated (see Fig. 11 in Stricklej 
f 19981 ) and the movie at http://www.uwm.edu/~jrs/out_of_corners.htm). 
Thus the motion mimics that of a large number of animals in a central at- 
tractive potential, but in this case arises solely due to a long range, indirect 
attractive interaction between the animals. Zooplankton in the dark, or in the 
aforementioned aquarium without a single individual in the beam to scatter 
light, show no evidence of a direct long range interaction. 

In the theory we can add the attractive interaction summed over all particles 
and this will enter an equation of motion similar to Eq. (1) but now repre- 
senting the balance of forces on one of the particles, i, in the ensemble, 

md t Vi = -70 Vi + d e^t) - K h 



and the energy uptake of the i-ih particle is, 

d t ei (t) = q - cei (t) - d v\ d (t) . (8) 



r,- 



V 



(7) 



Equations (3) and (4) remain the same except that the velocity must be in- 
dexed: Moreover, if we build in the effect of population variability, then the pa- 
rameter, d — > di in Eqs. (3)- (5), where the population noise in Eq. (6) results 
in a Gaussian distribution of the di. Notice that in Eq. (7) the term in square 
brackets has replaced the force generated by external potential — VU if). More- 
over, the force is linear in the distance between particles, r*j — fj, so that the 
mean field potential is quadratic as before, and its strength is adjusted with 
the coupling constant k^. 

The many-particle theory as outlined here, leads to results similar to those 
predicted by Eqs. (3)- (5), in particular, there is a bifurcation from a noisy 
fixed point motion to the symmetric pair of limit cycles except now the entire 
population participates in these motions. For simplicity, we treat this problem 
only considering the environmental noise, appearing in Eq. (1), though 
in the many particle case, Eq. (7), an independent noise, applies to each 
particle in the population. The appropriate measures are thus probability den- 
sities in both coordinate and velocity space. In addition to these motions, the 
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theory predicts the third motion of our suite: swarming. A large population 
of particles, initially widely and randomly dispersed, as shown in Fig. 11(a), 
begin to migrate toward the center of mass of the colony when the inter- 
particle interactions are switched on as shown in Fig. 11(b). Their individual 
motions are largely rotational, but this is not visible to the eye, since on aver- 
age equal numbers rotate in opposite directions. The CW and CCW motions 
are marked by arrows in the Figure. Figure 11 shows example trajectories but 
the appropriate measures are probability distributions of spatial and velocity 
coordinates on the x-y plane. 




(a) 



(b) 
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Fig. 11. Swarming in the many particle ABP theory, (a) Initially many particles 
are widely dispersed to random positions, (b) Under the influence of their mean 
field, they move toward the center of mass, (c) Directed motion, (d) Noise induced 
collapse of translational velocity to zero mean translation for D > D cr it- (e) Both 
sets of initial conditions finally converge on an equilibrium state is established with 
approximately equal numbers of particles rotating CW and CCW. The mean radius 
of the two limit cycle distributions depends on the mean self-propelling velocity. 



The simulations of the many particle theory show another interesting phe- 
nomenon that is observable for a different set of initial conditions as shown in 
Fig. 11(c). The particles are still initially randomly distributed over the plane, 
but now each is given the self-propelling velocity and all vectors are initially 
parallel (in this case pointed in the x direction). We call this directed motion. 



16 



For noise intensity smaller than some critical value, D < D cr i t , the particles 
continue to move with a mean translational velocity, but for D > D crit , the 
mean translational velocity rapidly descends to zero and the motion is all rota- 
tional (again with equal probability for C W and CCW motions). This is a new 
type of noise induced transition (see e.g. iHorsthemke and Lefeverl . [1984) and 
it is shown in Fig. 11(d). The critical noise intensity, -D cr it, depends on the pa- 
rameters of the velocity dependent diss ipation, 'y(v), and has b een calculated 



in the limit of large population density ( Erdmann et al. . 2002b| ). For both sets 



of initial conditions, the final state is the set of two symmetric counter-rotating 
limit cycles, as shown by the example trajectories in Fig. 11(e). 



6 Experiments with large Daphnia density 



Swarms are often formed for different purposes, such as enhanced feeding, 
mating and offspring rearing, by self-propelled animals living and moving in 
three-dimensional environments. Besides this, many animals are thought to 
form swarms to confuse and av oid pred ators dHumphries and Driver 



Hamilton!. Il97ll IPulliaml. Il97.4 see also iParrish and Edelstein-Keshetl . _ 
and lOkubo and Levin . 1200 if ) . Swarming for this purpose seems common to a 




number of anim als including planktivore fish, for example, sard ines, mackerel 
and anchovetta (|Partridgel.ll982l : lHall et all Il986f ) , zoopla nkton djakobsen et al. 
1994 ; Kvam and Kleivenl . ll995| ) and some species of birds ( Caraco et al.Lll980h . 



Such swarms, occasionally also associated with rot ational motions, have been 
called self-organized or emergent structures (see e.g 



1999HCamazine et alllioOlllOkubo and LevirlEoOlf) meaning that there is no 



'arrish and Edelstein-Keshet 



single leader nor a dynamical process that governs the motions of individual 
animals in such a way as to give rise to the swarm structure. 

We can consider swarming to be simply a local increase in density of animals 
without a net global c oherent motion. An attractant , such as a light mark er is 
usually necessary (e.g. Jakobsen and Johnsen . 19881 : Buskev et al. . 19961 ) . An 
additional collective motion, however, has been observed in swarms. According 
to anecdotal evidence, transitions from more-or-less random motions to an av- 
erage, coherent rotational motion have sometimes been observed in the field. 
These are called vortex motions a nd have actually been de scribed in detail 



for a species of oceanic plankton ( Lobel and Randall . 19861). In the labora- 



tory, yortex motions have been studied in b acterial colonies dBen- Jacob et aL 



11994 ICzirok et all Il996h and slime molds (|Rappel et all Il999f ) but have not 



been well investigated for the larger animals largely due to limitations imposed 
by the size of the swarms, for example of birds and fish. The physical, bio- 
logical, and chemical reasons for vortex-swarming are, consequently, not well 
understood. But the theories discussed here can provide at least a minimal 
set of conditions necessary for vortex formation. The symmetry of the pair 
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of limit cycles must be broken in order for the vortex to form. Symmetry- 
breaking processes based on pair avoidance (arising for example from short 
range repulsive potentials) are discussed in Sect. 7 below. In Sect. 8, we also 
consider hydrodynamic feedback as a symmetry breaking mechanism. 

We have observed transitions to vortex motions in Daphnia swarms using the 
apparatus shown in Fig. 7 as well as in a cylindrical aquarium (about 50cm 
diameter) recorded with a digital camera from above. The vertical visible light 
shaft, provided either by a flashlight beam or by the cylinder of Lucite, is the 
central attractant. First, Daphnia are conditioned in the dark for typically 
15 minutes in order that they assume an approximately random distribution. 
Then the light shaft is switched on. The following motions are typically ob- 
served. During the first few seconds, the animals swarm toward the light shaft 
as shown in Fig. 12(a). During this period no rotational motions can be ob- 
served; the animals simply rush toward the attractant. In the next period 
rotational motion becomes evident but with approximately equal numbers of 
individuals traveling in both directions (see Fig. 12(b)). Finally the individuals 
tend to align their velocities, and global rotational motion either in the CW 
or CCW direction occurs as shown in Fig. 12(c). We have observed that the 




Fig. 12. Three top view snapshots taken at 5 second intervals show the formation of 
a Daphnia vortex-swarm, (a) Daphnia gathering at the light shaft, taken 5 seconds 
after switching on the visible light, (b) rotational motion of swarming Daphnia in 
both directions, CCW and CW, and (c) vortex-swarming in CW direction. 

times for these processes to take place strongly depend upon the Daphnia den- 
sity, being shorter for larger densities. Using some inert particles to decorate 
the water motion, we have noted that not only do the animals move but also 
the water moves in the same direction after the vortex is fully formed. Thus 
we call this a biological- hydrodynamic vortex. Below in Sect. 8, we describe 
this phenomenon as a physical phase transition. Within the framework of the 
RWT, we hypothesize that the symmetry of the limit cycle pair is broken by 
a positive feedback represented by a non-zero order parameter that evolves 
in time. If a momentary fluctuation results in even a small majority of ani- 
mals circling in one direction with water being dragged along with them, then 
that direction becomes dominant as the other animals feel the hydrodynamic 
flow. Animals initially moving counter to the dominant flow then change their 
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direction thus forming the vortex. 

In Section 7 below, we show that in order for the vortex to form, some symme- 
try breaking mechanism must be present. We suggest three possible processes: 
short range repulsion, avoidance and velocity aligning hydrodynamic coupling. 
Daphnia do exhibit avoidance events as discussed below in Sect. 7. It is impor- 
tant to understand that avoidance is different from collisions between particles 
mediated by a symmetric short-range repulsive potential. 



7 How is the symmetry broken? 

As discussed above, the two motions characteristic of small density experi- 
ments with Daphnia, and described by the theory for non-interacting particles, 
are the noisy fixed point followed by a bifurcation to a pair of closed paths 
(limit cycles) that describe the equally probable CCW and CW rotational 
motions in a central potential. Many particles with attractive interactions 
can exhibit the third motion - swarming - this is represented by collective 
motions toward a center together with the aforementioned bi-directional rota- 
tions. The forth motion - the vortex - is also a many-particle phenomenon but 
is characterized by the single rotational motion of all the particles, and in the 
experiment, includes also fluid rotation. As we have seen, all these motions are 
observable in our apparatus, the former two with low Daphnia density and the 
latter two only in experiments with large densities. We also have mentioned 
that theoretically some symmetry breaking process is necessary in order to 
replace the two counter rotating limit cycles with a single vortex rotating in 
a single direction. The question is, in natural phenomena involving swimming 
animals, what is the symmetry breaking force or process? 

We put fourth two possibilities: short range repulsion and hydrodynamic self- 
alignment of velocities. Perhaps the simplest theoretical symmetry-breaking 
mechanism is offered by a symmetric short-range repulsive potential between 
pairs of particles. The physical laws governing scattering thus describe en- 
counters among such particles. Avoidance is, however, not the same as particle 
scattering. It involves animate behavior and arises, for example, in dynamical 
descriptions of pedestrian avoidance maneuvers. Hydrodynamic coupling in- 
volves the overlap of the flow field due to the motion of one moving particle 
with that of another. As we show below, this coupling tends to align velocities 
leading to a final rotational state in one direction. 

A third mechanism for symmetry breaking is collision avoidance. The first 
and second are discussed below in some detail, while this is only outlined as it 
has been described elsewhere. As we use the term here, avoidance has behav- 
ioral implications and thus involves mechanisms that cannot be described by 
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simple, symmetric potentials ( Mach et al. . 2003f ) as, for example, the Morse 
potential used below. We do not discuss avoidance in detail here. A similar 
problem arises in dyn amical descriptions of pedestrian avoidance maneuvers 
(jHelbing et all I1994T ) . These again lead to a symmetry breaking and conse- 
quently to various collective motions including the vortex formations of inter- 
est here. 



Both scattering due to a symmetric short-range repulsive potential and hy- 
drodynamic velocity alignment work equally well within the framework of the 
ABP theory to break the limit cycle symmetry. Unfortunately, at present the 
experiment does not provide enough information to distinguish between them. 
Moreover, Fig. 13, showing an actual encounter of two Daphnia, suggests that 
they execute a collision avoidance maneuver. Here two Daphnia are approach- 




Fig. 13. Track of a bottom view of an example avoidance event. The two Daphnia 
are shown in their initial positions, being at the same height in the water column (as 
only a 6mm thick horizontal layer is illuminated by an IR laser and line generator) . 
The tracking software is used to depict their trajectories through subsequent frames. 
In this example the Daphnia avoid each other, in a few other (rare) examples they 
actually collide. 

ing each other in the horizontal plane, avoid each other and depart on different 
trajectories. Though we have observed numerous encounters between pairs of 
Daphnia, unfortunately, the detailed statistics necessary for realistic modeling 
have not yet been compiled. 



7.1 Short-range repulsion 



Let us assume that the interaction between the particles shows long-range 
attractive together with short-range repulsive behavior. As suggested by simi- 
larities to interactions between molecules, we use the Morse-potential to model 
this type of interaction. The resulting collective effec t is the formation of clus- 
ters, a well known effect f rom s olid state physics ([Morse and Stueckelbere . 
Il929t iMorsel . Il92flt iKittell . 1995). In the limit of large inter-particle separa- 
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tion, the Morse potential becomes approximately quadratic, that is, similar 
to the external and mean field potentials used in Sections 2 and 5. As those 
sections show, the self-propelling velocity together with a quadratic poten- 
tial lead to rotational motions (the two symmetric limit cycles). But now, 
the short-range parts of the Morse potential parallelizes the neighbors and 
breaks the symmetry of the limit cycles, leading to rotation of the entire clus- 
ter in one (or the other) direction. As we know from the previous sections 
rotations are stable within such potentials. Thus we can think of a cluster 
with an initial condition such that all stationary velocities of the particles are 
equal in magnitude and additionally all have the same direction. Increasing 
the noise within the system leads to a phase transition as shown in Fig. 14, 
where, for D greater than a critical value, -D cr it, the cluster stops its directed 
motion. The net motion of the cluster vanishes. Finally, this can lead to a 
rotation of the cluster. In this new type of noise-induced phase transition the 
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Fig. 14. Noise induced transition from directed to motion with zero mean trans- 
lational velocity but non zero mean angular momentum (short range repulsive po- 
tential - Morse potential). 

mean translational velocity crashes to zero for D — > D crit and is replaced by 
a mean rotationa l velocity. The s t ability of the rotational motions was dis- 
cussed in detail i nlErdmann et aL ( 2002a ) and the noise- induced transition in 
Frdmauu et all (j2002hh. 



7.2 Hydrodynamics in the ABP theory 



Daphnia live in a low Reynolds number fluid environment. Hence, we can imag- 
ine that laminar flows of water mediate the interaction between individuals. 
The approach to use hydrodynamic interacti on terms to model this coupling 
has a l ready been used be fore for some species (|Weihsj .l ll973atlHuth and Wissell . 
1990; iNiwai Il994j , I l996h . With this in mind, we consider the dipole-like flow 



around one particle, j, that is moving at the self-propelling velocity. We want 



21 



to know the flow field of j at a distance at the location of particle i, as 
shown in Fig. 15. Considering similar contributions from all the other particles, 
a net hydrodynamic velocity Vp is present at the location i. 



A A 




Fig. 15. Geometry of the hydrodynamic coupling showing parallel and perpendicular 
velocity components of the animals, before coupling through the flow field (left) and 
after coupling (right). F = Kpvp, where the vector velocities are defined in Eqs. (11) 



In order to add a hydrodynamic interaction into Eq. (1), this stochastic equa- 
tion of motion can be rewritten in the following way: 



md t Vi = -livi) Vi + k f v f - VZ7(fj 



(9) 



where the term up vp is generated by all the other moving particles, and the 
parameter, Kp measures the strength of the hydrodynamic interaction. This 
term tends to align the direction of the velocity v with that of Vp. Thus particle 
i is finally dragged in the direction of motion of j. The aligning force, given 
by the new term in Eq. (9), originates with the viscosity of the surround- 
ing liquid, and ca n be modeled by a well-kno wn term in the Navier-Stokes 
equation (see e.g. Landau and Lifschit A 1987 ). Within the laminar regime, 
and for low Reynolds numb er, the flow can be modeled by additive Oseen- 
contributions (Oseen, Il915affl . This way of modeling hy drodynamic interac- 
tions is similar to the one used in the theory o f electrolytes dFalkenhagenl. 1971 ) 
and macromolecules faubbard and Douglasl 19931: Timoshenko and Dawson . 
199.4 iNeelov et all \200A Within this approx mation the flow field becomes 



R 



5 + 



valid for ^> R 



(10) 



where R is an effective hydrodynamic radius of the particles, and ® is a 
tensor product. Strictly speaking the Oseen-law is valid only asymptotically for 
Tij ^> R. We do not claim that this is the only way to introduce hydrodynamic 
interactions, nevertheless the Oseen-law has the right symmetry properties. 
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This is to be seen from the alternative way of writing Eq. (10): 
R . /? ( r 



VF(ri 



E 



' 13 ' v j j 



valid for 3> R . 



This shows that the Oseen-flow consists of parallel and radial components. 
The first term on the right of Eq. (11) is the component of the flow field 
parallel to velocity of the particle i, and it is this component that tends to 
align velocities of % and j. The second term is the radial component and tends 
to act along the direction of the line of sight between i and j. The alignment 
effect that results in symmetry breaking and thus large-scale coherent motions 
is due to the parallel components. 

If one imagines that we start with a population of active particles in a liquid 
medium, confined by a quadratic potential, due to the combination of active 
friction and external confinement on average half of the particles will rotate 
CW and the other half CCW. In contrast to particles that move in a vac- 
uum, coupling is mediated through the liquid. This leads to an interaction 
as described, for example, by Eq. (11). Because the hydrodynamic interaction 
aligns the velocity vectors of the individuals the dual rotational motion breaks 
down and the final motion is of the entire population rotating in one direc- 
tion. This direction is determined by random initial conditions. Additional 
motions, including clustering, global rotation of clusters and final rotational 
direction selection, res ult from this hydrodynamic in teraction. These have been 
described in detail in ( Erdmann and Ebelins . 20031 ). 



8 Random Walk Theory (RWT) 



8. 1 Introduction 



Discrete Random Walk Theories (RWTs) have been used extensively in biol- 
ogy, mainly to describe the irregular motions of bacteria (Berg, 1993f ). RWTs 
are iterative, that is a certain set of rules is applied to a particle located at 
a certain point in space in order to compute its future location and this pro- 
cess is repeated (iterated) to compute sets of jumps to successive locations. 
Since Daphnia naturally move in sequences of hops separated by pauses with 
a new direction of motion chosen at the end of each hop (see Fig. 6), RWTs 
can generate descriptive models of their motions. In this section we describe 
such a model and compare its predictions with experimental data. The ex- 
perimental results reported in Section 4 for Daphnia at low densities indicate 
that the dominant interaction is with the central attractant and that inter- 
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animal interactions are much less important for low animal densities. Con- 
sequently theoretical models of self-prope l led, interacting animals a nd their 
variants ( Vicsek et al. . 19951 : Levine et al. . 2001 : Couzin et al. . 2002) are not 
applicable to experiments at low densities with our animals. Though the ABP 
theory outlined in Section 2 adequately explains the noisy fixed point and 
limit cycle motions, we outline here a theory based on random walks. This 
theory has the advantage of being simple and easily simulated numerically, 
since it is only necessary to iterate and update an algorithm rather than to 
integrate a set of stochastic differential equations. 



Any theory of Daphnia motion at low densities within an attracting field must 
predict the two observed motions: noisy movements near and through a fixed 
point and the rotational motions about the attractant. An obvious question 
arises: what are the minimum ingredients necessary to describe these motions? 
To answer this question, and to simulate the observed behavior of single Daph- 
nia, we developed a self-propelled particle model based on random walks with 
the aim of being as simple as possible. The final model for low animal densities, 
closely related to the ABP theory, consists of only two ingredients: a short- 
range temporal correlation to simulate the general movement of the Daphn ia 
in darkness and an attraction to the light marker ( Ordemann et al. . 2003bl ). 



8.2 The Turning Angle Distribution 



Discrete random walk models having memory, or correlation, are powerful 
tools to simulate various biological, chemic al and physica l processes that are 
more complex than pure Brownian motion (Hughes], 19951 ). In the case of our 
Daphnia, a short-range correlation is indicated by the fact that the animal's 
choice of turning angle, a, at the end of each hop is not completely random 
but instead is correlated to some degree with the previous direction. We have 
observed a large number of hops of Daphnia swimming in the dark. From these 
data we have assembled the distribution P(a) of turning angles (DTA), which 
turns out to be a bimo dal symmetric di s tributi on with a pair of maxima 
at approximately ±35° ( Ordemann et all l2003bf ). Shown in Fig. 16 is the 
distribution -P(|a|) of the absolute turning angles, whic h resembles the results 
found for the DTA of the oceanic zooplankton copepod ([Schmitt and Seuront , 
2001). The presence of these maxima signals the existence of a short-range 
temporal correlation. If the directions of motion were completely uncorrelated, 
as for purely Brownian walkers, the DTA would be uniform. Thus Daphnia 
prefer moving in the forward direction and the most probable turning angles 
a max = ±35° may have evolved to optimiz e the amount of territory cove red 
while foraging in patchy food distributions ((Gerritson and Stricklerl . 119771 ). 



In constructing our RWT we introduce the random process (noise) by ran- 
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domly choosing a new turning angle, a, from the experimentally measured 
DTA shown in Fig. 16 for every iteration, or hop. This represents a temporal 
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Fig. 16. The experimentally measured distribution of absolute turning angles, \a\. 
The histogram was assembled from analysis of 1599 hops taken from eight animals 
moving freely in a dark aquarium with uniformly distributed food. The Daphnia 
were conditioned in the dark for a minimum of 15 min. before recording data. 
The turning angles were projected ont o the two-dimensional horizontal plane. Plot 
adapted from Qrdemann et al. ( 2003bh . 



correlation of only one iteration (time step). Though Daphnia may have a 
memory extending over a few hops, our observations indicate that this does 
not dominate the motion as the first hop memory does. Moreover, including 
only a one step memory keeps the model simple, which is our major aim. 
The particle in our RWT hops a fixed distance, Ax, in a time step, At, upon 
each iteration. Thus the magnitude of its velocity is a constant, and this is 
the self-propelling velocity. These model conditions are realistic, since the hop 
length and mean vel ocity of Daphnia are observed to be approximately con- 
stant (|DodsorJ . E996Tl . 



8.3 Motions Stimulated by Light 



In order to simulate the behavior of Daphnia in a light field an attraction 
to the marker must be built into the theory. The ABP theory detailed above 
models the attraction with a linear restoring force in the direction of the center 
of motion. Similarly, we construct a "kick" proportional to the distance, r, of 
the particle to the light source and proportional to a normalized strength 
parameter L/{1 — L), where [0 < L < 1]. For every time step, the walker 
randomly chooses a directional change, a, from the DTA, then the kick, [L/ (1 — 
L)]r towards the light is applied resulting in a new heading angle 9. The final 
heading vector, labeled t i+ i — ti, is rescaled to unit length to maintain constant 
velocity of the walker. The geometry of a pair of iterations for the time steps 
At = ti — and At = t i+ i — ti is shown in Fig. 17. 
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Fig. 17. The geometry of a pair of iterations, ij and ij+i, of the RWT showing 
distance to the light source r, heading angle 6, and turning angle a, as taken from 
the DTA shown in Fig. 16. 

With the RWT outlined above, the moves of the walker can be simulated for 
varying attraction strength L. For L = 0.4 the walker revolves around the 
light source in both directions and frequently changes its rotational direction 
similar to the recorded Daphnia motions shown in Fig. 8. Thus the RWT 
predicts limit cycle motions similar to those of the ABP theory. For smaller 
values of L, the circular motion is less pronounced as the influence of the 
randomness increases. For larger L, the circling breaks down and the walker 
mainly steps back and forth over the light source. This simulates the noisy 
fixed point motions similar to those predicted by the ABP theory. But in the 
RWT, the strength of attraction to light governs the type of movement. That 
is, in the RWT, the shape of the potential varies while the self-propelling 
velocity is constant. Whereas, in the ABP theory, the shape of the potential 
is fixed while the self-propelling velocity varies (this velocity is determined by 
the parameters of the energy depot and the friction coefficient). 

We turn now to measures of average quantities as predicted by the RWT. As 
in the experiment discussed in Section 4, rotational motion is demonstrated 
by the distribution of heading angles, P{9). The model predicts a similar 
distribution with equally probable rotational motions in both directions as 
shown in Fig. 18(a) for several values of the light strength, L. The distribution 
of radii, P(r), is shown in Fig. 18(b) again for the same set of values of L. 
Fig. 18(c) displays the average number of iterations (Mc) before a change in 
the circling direction for different L. A maximum at L — 0.4 implies that 
the most pronounced circling takes place at this intermediate value of the 
attraction strength. Figure 18(d) shows the angular momentum distribution, 
P(L ang ) for this optimal value of L. Finally, the distribution P(Mq) of the 
number of hops (iterations) before a reversal of direction is shown in Fig. 18(e), 
again for L = 0.4. For comparison of the simulation results to the experimental 
data see the fig ures in Section 4. A mor e detailed discussion of these results 
can be found in Ordemann et al. ( 2003bj ). 
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Fig. 18. Statistical measures predicted by the RWT for L = 0.0 (black squares), 0.2 
(dotted line), 0.4 (black circles), 0.6 (dashed line), and 0.9 (black diamonds), see 
also inset of (b). (a) The distribution of heading angles, P(6). The twin maxima at 
±90° demonstrate equally probable rotational motions in both directions (CW and 
CCW). (b) The distribution of radii, P(r). (c) The average number of iterations 
{Mq) before a change in the circling direction takes place for different L. Note 
the maximum at L = 0.4 indicating the case of most pronounced circling, (d) The 
distribution of angular momentum, P(L ang ), and (e) the distribution of the number 
of hops before a change of rotational d irection, P{Mq), both at L = 0.4. Plots are 
adapted from lOrdemann et al. | (|2003hl l^. 
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8.4 Diffusion in RWT I 



We can study the natural diffusion of Daphnia by once again considering their 
motions in the dark. Within the framework of the RWT, a convenient measure 
of diffusion with time, is the end-to-end distance R. This is the straight-line 
(Euclidean) distance from the starting point of a walkers trajectory to the end 
point at the time £j. Of course, an ensemble average of such distances must be 
obtained from a large number of trajectories, each generated with a different 
seed for the random number generator. Of interest is the exponent, u, defined 
by (R 2 ) 2 oc t u , where t is the time (or the number of iterations) after the walker 
is started, and the ang ular brackets i ndicate the aforementioned ensemble 
average (see for example iHug hel Il995h . For purely Brownian particles, that 
is those with zero correlation, v = 1/2. In contrast, for our particles with a 
temporal short-range positive correlation, we find v = 1 for short time scales, 
followed by a crossover to v = 1/2 for asymptotically long times. This is 
shown by the simulation results (full circles) in Fig. 19. In general, for short- 




Fig. 19. Diffusion of Daphnia and of random walkers with and without short time 

i 

correlations. The root mean square (rms) end-to-end distance {R 2 ) 2 is plotted 
versus time, t. For the simulations, 500 realizations were used with a one step corre- 
lation (full circles) and zero correlation (open circles). The experimental data (open 
squares) was obtained from the trajectories of eight Daphnia, where the distances 
are in mm and the time is in sec. The slope of the dashed line is 1 and of the solid 
lines 1/2, indicating ballistic (correlated) and uncorrelated motions, respectively. 

range temporal correlations, the crossover time from non-Brownian motion to 
Brownian motion depends on the magnitude of the correlation time (memory). 
Here the correlation time is just one time step, since we simply update the 
turning angle according to the DTA at each iteration. For comparison, we show 
a simulation for purely Brownian walkers (open circles in Fig. 19) for which v = 
1/2 always. When calculating the rms end-to-end distance of Daphnia motions 
in the dark, we observe an exponent v = 1 for t < 30sec (see squares in Fig. 19). 
Thus their motions are gov erned by positive correla tions at least for short time 



scales (for more details see lOrdemann et all 12003(1 ). It is widely assumed that 
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for long time scales, that require much larger water reservoirs, the tracks 
of zooplankton in the field resemble Brownian motion, in agreement with 
our simulation at long times. Unfortunately, we were unable to observe the 
crossover experimentally, because we cannot record long enough trajectories. 
The trajectory length was limited by the size of our aquarium and by the 
available field of view in respect to the resolution of the video cameras. 



8.5 Diffusion in RWT II 



As has been observed in the experiments Daphnia tend to turn with a preferred 
angle at every hop they make (see Fig. 16). The question which occurs, having 
this in mind, is: Why do Daphnia hop with an DTA as shown above? What 
could be an reason for it? Can it be an evolutionary one? In order to answer 
this question it is reasonable to calculate the diffusion coefficient, having the 
aforementioned DTA i n mind, and c o mpar e this with normal Brownian motion 
as it has been done in Komin et al. (120041) . In order to calculate the diffusion 



coefficient we derive the mean square displacement of a swarm after n = t/ At 
time steps: 

n—l n 

(R 2 n ) = nAx 2 + 2Ax 2 ]T ]T (cos a) , (12) 

8=1 j>l 

which depends on the angular correlation function: 

7T 

C a = (cos a) = J P(a) cos a da . (13) 



After some transformations ( Okubo and Levinl . 2001 : Wu et al. . 2000j ) of Eq. (12) 
the mean squared displacement of a two-dimensional cor related random wal k 



can be calculated what leads to the diffusion coefficient ( Komin et al. . 2004 ): 



AD 



1 + C a Ax 2 
l-C„ ~At 



(14) 



As Eqs. (12) and (14) show the angular correlation function is crucial for 
the behavior of the time evolution of the mean square displacement of the 
particles and therefore for the diffusion behavior. If < C a < 1 the diffusion 
would be enhanced. This could be compare d with a mean turn ing angle < 
a < 7r/2, as observed with the Daphnia. In iKomin et al.l (J2004J) C a has been 
calculated analytically assuming the DTA can be approximated by a Gaussian 
distribution (shown in Fig. 20). For (a) = n/2 we have a C a = and therefore 
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(a) Angular correlation 



(b) Diffusion coefficient 



Fig. 20. (a) Angular correlation versus mean (a) and variance a a of a Gaussian like 
DTA. (b) Diffusion coefficient. The parameters out of the Daphnia experiments are 
indicated with an arrow at the following values: {a) = 0.84, a a = 0.63 and D = 3.1. 

D — 1, what is equivalent to free diffusion, and does not depend on the 
standard deviation of the turning angle anymore. 

Taking the experimental parameters for the Daphnia looking at Fig. 20(b) it 
seems to be obvious that the DTA for zooplankton has evolutionary advan- 
tages. The diffusion is enhanced and therefore an animal is able to cover a 
larger area while foraging for food. 



8.6 Breaking the Symmetry in the RWT 

We now turn to a simulation of the transition to vortex viewed as a phase 
transition within the context of our RWT model. The model is modified based 
on ideas about hydrodynamic feedback. We assume that the Daphnia are 
sensitive to small water movements, and that they have a tendency to align 
the direction of their motions with that of the majority of their neighbors. We 
thus define a global order parameter, O = (New — -^ccwV-^tot , where A^cw 
and A" ccw are the instantaneous number of particles moving CW and CCW, 
respectively, and N tot is the total number of particles. Particles are considered 
to be executing rotational motion if the absolute value of their heading angles 
lies within the range, 80 < \9\ < 100. In the simulation, O evolves with 
time. It is zero at instants when there is no net directional motion of the 
particles, smaller or larger than zero when there is net CCW or CW motion 
respectively. For the fully formed vortex (all particles moving in the same 
direction), O — > ±1. 

We now focus on motions in a small neighborhood of radius a, supposing that a 
single Daphnia can sense hydrodynamic disturbances only from its neighbors. 
Motions within the neighborhood of a particle are governed by a local order 
parameter, V = (Aqw — Accw)/N a , where Aqcw and Aqw are the numbers 



30 



of particles within the neighborhood that satisfy the same definition of rota- 
tional motion, and N a is the number of all particles within the neighborhood. 
In addition to the kick toward the light that applies to all particles in the 
simulation, we now include a local kick, K = qV in the appropriate direction 
(CW for V > 0, CCW for V < 0) to the corresponding particle. Depending on 
the strength, q, of the kick and the size of the neighborhood, a vortex forms 
after a certain time. The formation of the vortex can be tracked by observing 
the temporal evolution of the global order parameter O. An example is shown 
in Fig. 21. 




50 150 25C 50 150 250 



t t 

(a) (b) 

Fig. 21. Evolution of the global order parameter, O, showing the formation of 
the vortex (a) in CW direction (O — * +1) at t ~ 190 time steps for parameters, 
N tot = 30, a = 0.7, and q = 0.16, (b) in CCW direction (O -» -1) at t ~ 180 time 
steps for the same parameters. 



9 Discussion and Conclusions 



Daphnia provide for the first time, the opportunity to construct well-controlled 
laboratory experiments with real biological agents. Such experiments provide 
crucial tests of existing different theoretical models of swarming, fixed point 
and rotational motions and vortex formation. And therefore these experiments 
open up the possibility to learn more about the aforementioned four differ- 
ent animal and agent motions through the interplay of theoretical modeling 
and experimental investigations. For Daphnia, the interaction that leads to 
swarming and vortex formation is not a direct one, as is the case for birds and 
fish that almost certainly use visual velocity alignment. Moreover, Daphnia 
swarms do not perform ob yious predator ay oidance maneuvers such as swarm 
splitting and recombining ( Hall et al. . 1986f ). Nevertheless, the observation of 
vortex-swarming by Daphnia in the laboratory provides the possibility to learn 
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more about the physical, chemical and biological aspects of this phenomenon 
that may be important for understanding similar motions exhibited by other 
living creatures. Further experiments with Daphnia will be necessary to un- 
derstand in detail the factors that influence vortex-swarming, for example, 
Daphnia species, predator kairomone concentration, density of both food and 
individuals, and light intensity and wavelength. Light patterns and their per- 
ception by Daphnia as well as the physical aspects of the bio-hydrodynamic 
vortex need more attention. 

The ABP model assumes that particles can absorb energy distributed over a 
medium in which they move and store it in a depot in a way similar to animals 
actively foraging for food over some space. Like animals, the particles must 
expend energy in order to move through the medium, and they experience a 
velocity-dependent dissipation. Also they expend energy to stay alive (similar 
to metabolism). Furthermore, also like real animals, they are subject to vari- 
able (random) forces, or "noise" . A confining potential is added in order that 
the particles are confined (or swarm). The model is therefore well motivated 
from a biological point of view. Interestingly, in addition to swarming, it pre- 
dicts two basic motions that we designate the noisy fixed point and the equally 
probable, symmetric limit cycle pair. Similar motions are commonly observed 
in a variety of animals ranging in size and complexity from bacteria to birds 
and fish. Moreover, the ABP model has taught us that these two canonical 
particle motions arise from a minimal set of conditions. The particles must 
exhibit a mean self-propelling velocity and they must be confined either by an 
external potential or by the mean field that arises from an interparticle attrac- 
tive interaction. Our experiments with small densities of Daphnia demonstrate 
the two basic motions. Moreover, the two limit cycles form the skeleton from 
which the vortex is built after some symmetry-breaking process eliminates one 
in favor of the other. At large particle densities the APB model was extended 
to include symmetry breaking by Oseen hydrodynamic flow in one instance 
and by an interparticle avoidance process, modeled by including a short range 
repulsive potential, in another. Both processes lead to vortex formations, sim- 
ilar to those observed with high density Daphnia swarms, but it remains to 
investigate which one better describes the experimental observations. 

The model described by the RWT has also taught us something significant. 
Without including the peaked DTA, that signifies a statistical bias in favor of 
forward motion, in the model, no limit cycle motions of the walkers develop 
for any value of light intensity. With a uniform distribution of turning angles 
(uncorrelated Brownian walkers), only noisy fixed point motion is possible. 
Thus short-range temporal correlations, here indicated by the twin peaks in 
the DTA, are necessary for rotational motions about the central attractant. 
Within the framework of the RWT, we can now list four minimal conditions 
necessary to form vortices: i) a self-propelling mean velocity, ii) a statistical 
bias in favor of forward motion, iii) a central attractant, and iv) a symmetry 
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breaking process that tends to align local velocities. The bimodal DTA, nec- 
essary for rotational motion, might be an essential ingredient of the collective 
motions of large colonies of creatures and thus might be a general feature 
detectable in other animals as well. Why the natural movement of plankton 
follows such a distribution with maxima located at approximately ±35° (true 
for both our Daphnia and ocean-going zoo plankton) and how the animals 
actually choose to turn either left or right at the end of each hop are unan- 
swered questions. But it is intriguing to speculate that these preferred turning 
angles may somehow enhance fitness, for example by possibly maximizing the 
acquisition of energy when food is distributed in patches. Indeed, we have 
calculated the diffusion coefficients for a swarm of random walkers using the 
experimentally determined bimodal DTA and found that, for parameters ap- 
propriate to the Daphniay, diffusion is considerably enhanced over what one 
would obtain for a walker with purely uncorrelated motions. This strongly 
suggests that motions characterized by a preferred turning angle confer an 
evolutionary advantage on foraging animals. 
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